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SUMMARY 


Populaticai parameters of 275 streams (includii^ 20 additional streams) in the 
synoptic-year sample have been foimd by a computer technique. Some 16% of the 
sample is in these streams. 

Four meteor streams that have close orbital resemblance to Adonis cannot be 
positively identified as meteors ejected by Adonis within the last 12000 years. 

Ceplecha's discrete levels of meteor height arts not evident in cur radar meteors. 

The spread of meteoroid frs^ments along their common trajectory has been com- 
puted for most of our observed meteors. There is an unexpected relationship between 
spread and velocity that perhaps conceals relationships between fragmentation and 
orbits; a theoretical treatment is necessary. 

Revised unbiased statistics of synoptic-year orbits are presented, together with 
parallel (new) statistics for 1961-65 orbits. 
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PHYSICAL AND DYNAMICAL STUDIES OF METEORS 
Meteor-Fragmentation and Stream-Distribution Studies 


Zdenek Sekanina and Richard B. Southworth 
Center for Astrophysics 

Harvard College Observatory and Smithsonian Astrophysical Observatory 
Cambridge, Massachusetts 02138 


1. DISTRIBUTION OF METEORS IN THE STREAMS DETECTED 
IN THE SYNOPTIC- YEAR SAMPLE 


1. 1 Introduction 


We reported previously (Southworth and Sekanina, 1973) that our two-phase com- 
puterized stream search among 19, 698 radio-meteor orbits of the synoptic-year sample 
resulted in the detection of 256 streams. 

The distribution of radio-meteor orbits in these 256 streams has now been studied 
in terms of the D-test, which measures the similarity of any two orbits by the differences 
in their Keplerian elements (Southworth and Hawkins, 1963). The D-test was utilized 
by Sekanina (1970a) in the formulation of his statistical model of meteor streams. The 
model describes the D-distribution function of meteors in a stream by two parameters; 
the population coefficient A, which measures the stream’s strength relative to the level 
of the ambient sporadic background; and the dispersion coefficient <J, which determines 
the degree of scatter in the distribution of meteor orbits within the stream. 

1.2 New Method of Determining the Distribution-Function Parameters A and (J 

In our previous studies of radio-meteor streams (Sekanina, 1970b, 1973), the two 
parameters of the D-distribution function were derived graphically from the log-log 
plots of the cumulative number of meteors versus the D-test (Sekanina, 1970a). To 
avoid the laborious plotting of hundreds of distribution curves, we have now developed 
a new, fully computerized method. 



The theoretical cumulative D-distribution of the statistical model of meteor streams 
has the form 



where 9t(D) is the number of meteors of the studied sample whose values from the 
D-test are less than or equal to D, relative to the mean orbit of the stream in question; 
C is the normalizing constant; and 


f(x) = x^‘ ^ + 2. 64 a[ r e ^ dt - xe ^ 


( 1 - 2 ) 


\0 / 

Since the parameters A and o cannot be determined explicitly, an iterative least- 
squares procedure has been proposed as follows. 


We start from an arbitrary pair of values for the two parameters, A^ and cr^, and 
calculate ?^q(D) from equation (1-1) with the proper C^. Then for each point on the 
D-distribution curve, we calculate tl^us establish a number of conditional 

equations of the type 


A log 7? H log 


JL= g-log K 

1\q 
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9 log (7 
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(1-3) 


where A log A = log A^ - log A^, . . . are the corrections to be applied to the param- 
eters to obtain their improved values A^, .... Substituting A^ for A^ and calculating 
new ^q(D) values from equation (1-1), the procedure can be iterated until it converges. 
The normal equations for the corrections to the parameters, implied by the conditional 
equations, are 
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where n is the number of points on the cumulative D-distribution curve, 

A(,= -2^f^i^(l+ 1.389 AE‘“-®e'®^ , 

A =1-^ 

^ f(E) ’ 


(1-5) 

( 1 - 6 ) 


and 


E = 
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(1-7) 


1. 3 The D-Distributions of the Detected Streams 


Since the cumulative D-distribution of meteors in a stream is a by-product of the 
stream-search program, the above differential-correction procedure can work directly 
with the punched list of meteors from the search program. 

The differential-correction procedure was first tested on several streams of the 
1961—65 sample, for which A and O’ had been previously determined graphically 
(Sekanina, 1970b, 1973). A few examples, listed in Table 1-1, indicate that the agree- 
ment between the least-squares solution and the graphical method is reasonably good. 

Practical calculations have shown that the selection of the initial values for A and 
o has no effect on the convergence of the differential-correction procedure. We chose 
A = 40 and 0 = 0. 050 as standard initial values and found that the procedure converged 
successfully for all but six of the 256 streams of the ssmoptic-year sample. The max- 
imum number of required iterations was 70, but in most cases less than 10 were 
necessary. In four cases (May Arietids, a Draconids, L Cepheids, and € Umids) the 
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Table 1-1. Parameters of the D-distribution from the least-squares 
and graphical methods. 


Stream 1961-65 


Least- 

-squares 

Graphical 


A 

a 

A 

a 

Quadrantids 

85 

± 19 

0.038 ± 0.002 

50 

0.037 

Lyrids 

154 

± 40 

0.039 ± 0.003 

200 

0.043 

October Draconids 

4.4 ±0.6 

0.064 ± 0.006 

5 

0.066 

Geminids 

189 

± 35 

0.037 ± 0.001 

200 

0.038 


procedure failed to converge, and in two (April Ursids and e Ursids) it failed to 3 deld 
any solution, because the slope of the D-distribution curve at D < 0. 1 was steeper than 
allowed by the model. In one other case (a Aurigids), the solution indicated that this 
"stream" does not satisfactorily discriminate from the sporadic background. The 
population and dispersion coefficients of the six streams for which the differential- 
correction procedure failed have been determined graphically. 

1.4 Numerical Results 


The results for the 255 streams (omitting a Aurigids) are listed in Table 1-2. The 
individual columns give the following information (calculated, in many instances, from 
the referenced formulas of Sekanina, 1970a): 


Column 1. Designation of the stream. 
Column 2. Population coefficient A. 


Column 3. Dispersion coefficient a. 

Column 4. Calculated cumulative number of meteors with D < 0. 20, 
according to Sekanina (1970a). 


0.20 


I 


N(D) dD 


Columns 5 and 6. Inner (or differential) Dj and outer (or integral) Djj limits of 
the stream [formulas (38)] . They correspond, respectively, to D-values of equal 
density and equal population between the stream and the background. 


4 


Column 7. Relative sporadic contamination (%) between D = 0 and 0. 20 [formula 
(33)] . 

Column 8. Stream-to-background concentration ratio at D = 0. 20 [formula (32)] . 

Columns 9-11. Calculated numbers of definite members of the stream within, 
respectively, D = 0.20, D = Dj, and D = Djj [formula (34)]. 

Column 12. Calculated total number of definite members of the stream [formula 
(35)] . 

Columns 13 and 14. p- and C-tests, respectively, of the existence of the stream 
[formulas (43) and (48)] . 

Colunrn 15. Cosmic weight of the stream, as defined by Whipple (1954). 

Colunms 16 and 17. The date (CST) of the middle of the period of activity and its 
duration (days). 

Columns 18—20. Observed cumulative numbers of meteors with D < 0.20, 

D < 0. 25, and D < 0. 30, respectively. The degree of consistency between columns 18 
and 4 is a measure of the differential-correction fit of the distribution curve. 
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Table 1-2. A list of D-distribution parameters of 255 streams of the synoptic-year sample. 
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2. ADDENDUM TO THE STREAM SEARCH AMONG THE RADIO METEORS 
OF THE SYNOPTIC-YEAR SAMPLE 



Since the distribution of meteors of the S 5 moptic-year sample showed a sharp drop 
at high inclinations and an extremely low level of population in retrograde orbits, we 
have relaxed somewhat, for these types of orbits, the restrictions in the screening 
procedure we used to delineate the 256 streams. 


The relaxed screening resulted in the detection of 20 more high-inclination and 
retrograde streams, which are listed in Table 2-1. This table is arranged in the same 
way as Table 9-2 of Southworth and Sekanina (1973). For the sake of convenience, we 
list here the contents of the individual columns (all angles are referred to the standard 
equinox 1950.0); 

Colunm 1. Designation of the stream. 

Column 2. Argument of perihelion w and its mean error (degrees). 

Column 3. Longitude of the ascending node n and its mean error (degrees). 

Column 4. Inclination i and its mean error (degrees). 

Column 5. Perihelion distance q and its mean error (a. u.). 

Column 6. Eccentricity e and its mean error. 

Column 7. Semimajor axis a (a. u.) and revolution period P (years). 

Column 8. Longitude of perihelion u = co + (degrees) and the date of passage 
through the node (UT and days in 1950). 

Column 9. No-atmosphere velocity V and its mean error (km sec ^). 

oo 

Column 10. Right ascension of the corrected mean radiant from individual 
meteor radiants and its mean error (degrees). 
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Column 11. Declination 6„ of the corrected mean radiant from individual meteor 

K 

radiants and its mean error (degrees). 

Column 12. Right ascension and sun-oriented celestial longitude - Xq of 
the corrected mean radiant from the stream's mean orbit (degrees). 

Column 13. Declination 5^^ and celestial latitude of the corrected mean radiant 
from the stream's mean orbit (degrees). 

Column 14. Geocentric V„ and heliocentric V„ velocity at the node (km sec ^). 

(jr ri 

Column 15. The node (A = ascending; D = descending) and the type of stream 
(C = circumpolar, with the mean radiant permanently above the horizon at Havana, 
Illinois; D = daytime; M = mixed, with the mean radiant never more than 10® below the 
horizon; N = nighttime). 

Column 16. Height at maximum ionization and its mean error (km). 

Column 17. Number of meteors in the sample with D < 0. 25 with respect to the 
mean orbit, with an asterisk to indicate that the stream was detected in the 1961-65 
sample. 

The additional search increased the number of retrograde streams from the 
previous two (Perseids and Orionids) to 15. The most abundant of the new retrograde 
streams, August Triangulids, is about as populous in our sample as the Orionids are. 
The other ones are more comparable in population to the Perseids, with an average of 
stream members of only three or four. 

2. 2 Revisions among the 256 Previously Reported Streams 

The considerations involving the distribution of meteors in the streams and their 
periods of activity resulted in several minor changes in the list of the detected streams. 
We already commented (Section 1. 3) that the analysis of the D -distribution of a Aurigids 
showed that the stream does not satisfactorily discriminate from the sporadic back- 
ground. 
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Table 2-1. A supplementary list of streams detected in the synoptic-year sample. 
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Table 2-1 (Cent.) 
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The periods of activity of five streams required adjustments that also brought 
about minor changes in the streams' mean orbital elements. The revised data on 
these streams are also included in Table 2-1, replacing the corresponding figures in 
Table 9-2 of Southworth and Sekanina (1973). 

2.3 The D-Distribution of the Additional and Revised Streams; Total Number of Stream 

Meteors in the Sample 

The data on the distribution of meteors in the five revised streams, as listed in 
Table 1-2, are based on the corrected figures. The distribution parameters of the 20 
additional high-inclination and retrograde streams are given in Table 2-2, which is 
arranged identically to Table 1-2. For the details, see Section 1.4. 

The minor high-inclination and especially retrograde streams have mostly rather 
low dispersion coefficients, typical for compact streams, and very high population 
coefficients, indicating practically no interference from the sporadic background in 
small values of D. 

In total, we have delineated 275 streams in the synoptic-year sample (256 originally 
detected minus a Aurigids plus 20 additional). The total number of definite stream 
members in the 275 streams is 3142, which is about 16% of the whole sample. The 
average population of the detected streams is 11 to 12 meteors per stream. 
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Table 2-2. A list of D-distribution parameters of 20 additional high-inclination and retrograde streams. 





3. THE POTENTIAL ASSOCIATION OF FOUR METEOR STREAMS WITH 

THE MINOR PLANET ADONIS 


3. 1 The Problem 


Four streams in the synoptic-year sample have mean orbits that match the orbit 
of the minor planet Adonis within D = 0. 20 (see Table 10-1 of Southworth and Sekanina, 
1973): X Sagittariids, Scorpiids-Sagittariids, « Aquarids, and Capricormds- 
Sagittariids. 

The orbital similarity may, of course, suggest an evolutionary relationship 
between the stream and the minor planet, which in turn could imply that Adonis might 
have been an active comet long ago and that the associated meteor streams currently 
observed represent its debris. In that case, the streams should indeed move in orbits 
similar to but not identical with that of Adonis. The orbital difference is partly due to 
the nonzero momentum gained by the meteoroids at ejection and enhanced later by the 
accumvilation of differential perturbations by the major planets, and partly due to 
the effects of radiation pressure on the meteoroids. The Pojmting-Robeitson effect 
can be shown to accumulate, over the spans of time considered here, to no more 
than 0.1 a.u. in the semimajor axis, which is less than the uncertainty in the semi- 
major axis of the mean orbits of the streams. 

3. 2 The Calculations 


We have attempted to explore the observed difference between the orbits of Adonis 
and the orbits of the potentially related streams to learn something about the time and 
circumstances of ejection. 

Since Adonis has its aphelion at 3.3 a. u., close encounters with Jupiter are 
excluded. Because of this circumstance and the rather low precision of the starting 
data, we considered only secular perturbations in order to facilitate calculations as 
much as possible. 
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Applying the same method that was used to investigate a potential evolutionary 
relationship between Adonis and a radio-meteor stream from the 1961—65 Smithsonian- 
Harvard sample (Sekanina, 1971), we started our model calculations by running the 
orbit of Adonis for 12000 years backward in time, taking the secular perturbations by 
Jupiter to Neptune into account. Next we varied the following quantities to study the 
effects on the orbital modification of a meteoroid ejected from Adonis: 

A. The time of ejection (on a gross scale). 

B. The position of Adonis in orbit at the time of ejection (relative to perihelion). 

C. The magnitude of the velocity of ejection. 

D. The direction of ejection (relative to the sunward direction in the orbit plane). 

We adopted five discrete times of ejection, 4000 to 12000 years ago, 2000 years 
apart. For each time, the ejections were assumed to take place at 1. 2 and 0. 7 a. u. 
from the sun before perihelion, at perihelion, and also at 0.7 and 1.2 a. u. after 
perihelion. Furthermore, five different directions of ejection were considered: toward 
the sun and deviating, both ahead of and behind the solar direction, by 30° and 60°. 
Finally, the magnitude of the ejection velocity was estimated from Probstein's (1968) 
fluid-dynamics model. Assuming that at the time of ejection the nuclear radius of 
Adonis was between 1 and 20 km and that the surface was covered by water snow (con- 
taminated by meteoric matter), and accepting that at 1 a. u. from the sun the vaporiza- 

17 17 -2 -1 

tion rate was between 1 X 10 and 7 X 10 molecules cm sec (owing to uncer- 
tainties in the emissivity of the surface), we can establish from Probstein's model 

-4 

that the ejection velocity of meteoroids of individual masses of about 10 g and 

-3 -1 -1 

densities of 1 g cm should vary between 17 and 200 r (m sec ), where r is the 

heliocentric distance in a. u. Calculations were carried out for the two extreme values 

as well as for 60 r ^ (m sec ), which has been thought to be a reasonable mean 

value. 

3. 3 Corrections to the Initial Elements of Ejected Meteoroids 

The magnitude of the impulse on a meteoroid separating from the parent body 
shows up in the differences between the orbital elements of the two bodies at the time 


24 


of ejection. If the known position of ejection in the orbit is given by the true anomaly 
V and the known velocity of ejection is defined by its radial component | (counted posi- 
tive away from the sun) and transverse component t) (positive in the direction of motion) 
in the orbit plane, we need to apply the following corrections to the orbital elements 
of the parent body: 



- 4 cos V + T) 


2 + e cos V 
1 + e cos V 


sin V 


> 


An = Ai = 0 , 


, . 2 
1 - cos V + e sin v 


Aq = -r^ 1-4 sinv + 7i ^ 
^l+e\^ ' 1 + e cos V 


Ae = B 


I sin V + T) 


• 2 cos V + e(l + cos^ v) 


1 + e cos V 


(3-1) 


where 

u 

_ 4 

and U = 2. 978 X 10 , if ^ and r\ are expressed in meters per second. 


(3-2) 


Furthermore, additional corrections have to be applied because of the effect of 

2 2 

radiation pressure. If ^ /k < 0 is the relative reduction in the gravitational constant 
due to radiation pressure, the corrections are 

2 

. _ _ sin V Ak 

rp " e ’ 


An = 

O 

II 


rp 

rp 


II 

1 - cos V 

Ak^ 

^ 1 + e 

k2 

Ae = 
rp 

- (e + cos v) 

Ak^ 
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Here v is the true anomaly at ejection and 


^ ^ 6 X 10"^ Q 

k2 " ■ Vs 


(3-4) 


where p and a are the density and radius of the meteoroid, and Q = 1 is its scatter- 
s s rp 

ing efficiency for radiation pressure. For meteoroid masses of the order of 10 g, 

we find pa =“ 0. 03 g cm ^ and, therefore, -0. 002. 

s s 

3. 4 The Results 


With five ejection times considered, each with 5 X 5 X 3 = 75 options in regard to 
the circumstances at ejection, we have a total of 375 individual model orbits per stream. 
For each of them, the initial orbital elements (at ejection) were computed from the 
orbital elements of Adonis at the time of ejection by adding the specified corrections 
for the effects of ejection velocity [eq. (3-1)] and radiation pressure [eq. (3-3)] . The 
secular perturbations by Jupiter to Neptune were then applied to calculate the present 
orbits of the 375 model ejections. These orbits have been compared with those of the 
four meteor streams potentially related to Adonis. 

The results are summarized in Tables 3-1 to 3-4 for Scorpiids-Sagittariids, 

X Sagittariids, Capricornids-Sagittariids, and e Aquarids, respectively. The tables 
list the values of the D-test (in units of 0.001) between the model orbits and the mean 
orbit of the stream. Values exceeding D = 0. 25 (i. e., >250 in the tables) have not 
been printed. The time of ejection and the magnitude of the velocity of ejection are the 
two basic parameters. The location in orbit at ejection is printed as a subparametei 
for each time of ejection, and the following abbreviations are used; 1 = before perihelion, 
heliocentric distance 1. 2 a. u. ; 2 = before perihelion, 0. 7 a. u. ; 3 = perihelion; 4 = after 
perihelion, 0. 7 a. u. ; 5 = after perihelion, 1.2 a. u. The direction of ejection is 
expressed in terms of the ejection angle, which is printed as a subparameter to the 
velocity of ejection. The ejection angle is counted from the sunward direction and in 
the orbit plane. Its negative values indicate the ejections "ahead of" Adonis, i.e., in the 
general direction of the orbital motion of the parent body; the positive values represent 
the ejections "behind, " i.e., opposite, the motion of Adonis. If we choose to associate 
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the angle of ejection with the spin orientation of the parent body, in analogy to Whipple's 
(1950) interpretation of the nongravitational effects on cometary nuclei, the positive 
ejection angles would imply direct rotation, and the negative angles, retrograde 
rotation. 

Inspection of Tables 3-1 to 3-4 suggests that the results do not allow us to make 
any clear-cut conclusions as to the existence of any genetic relationship between the 
minor planet and the four streams. The minima in the D-test distribution, which pick 
up the best models, are very shallow and in many instances barely show any improve- 
ment over the value of the D-test between the orbit of the stream and that of Adonis. 
Furthermore, the amplitude of the D-variations systematically increases with the 
velocity of ejection. Thus, we could formally improve the fit by increasing the ejection 
velocity to completely unrealistic values. Consequently, the ejection velocity cannot 
be determined from the variations in the D-test, but must be postulated. Therefore, 
in the tables we have marked the minimum values of the D-test separately for each of 
the adopted velocities of ejection. Even then, however, the age of the streams and 
the location of ejections in orbit are rather indeterminate, and only the direction of 
ejection comes out consistent for each stream. Scorpiids-Sagittariids, x Sagittariids, 
and Capricornids-Sagittariids show a preference for positive ejection angles, and 
e Aquarids, for negative. 
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Table 3-2. D-test comparison of model ejections from Adonis with the x Sagittariid stream. 
D(Adonis) = 0.089. (D's in units of 0.001.) 
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Table 3-3. D-test comparison of model ejections from Adonis with the Capricomid-Sagittariid stream. 
D(Adonis) = 0. 199. (D’s in units of 0. 001. ) 
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Table 3-4. D-test comparison of model ejections from Adonis with the € Aquaiid stream. 
D(Adonis) = 0. 193. p*s in units of 0. 001.) 
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4. HEIGHT-VELOCITY DIAGRAM 


4. 1 The Problem of Meteor Heights 

We reported previously (Southworth and Sekanina, 1973) that we had failed to detect 
in the synoptic-year sample the discrete levels of meteor height that had originally been 
reported by Ceplecha (1967, 1968) for photographic meteors. However, more recent 
inspection of the utilized height data indicated that the punched output of the main reduc- 
tion program (height-density cards), which served as the source, contained heights 
from diffusion for some meteors and geometric heights for others (plus radiant heights 
for a minor part of the sample) . The three types of heights differ from each other 
sometimes by as much as 10 km or more, and the geometric heights, which are gen- 
erally considered less reliable than the diffusion heights, constituted a significant part 
of the sample. To rectify this situation, we have now consulted the much more exten- 
sive printed (rather than punched) output of the main reduction program; we eliminated 
both the geometric and the radiant heights and have finally ended up with 14502 meteors 
with diffusion heights. 

4. 2 The Height— Velocity Plots 

Figure 4-1 shows one of the computer-generated plots of meteor height at maximum 
ionization h versus the no-atmosphere velocity Voo* We have also secured similar 
plots for beginning height; they look much the same, but have slightly larger scatter. 

Although the degree of precision of the present data is definitely greater than that 
of our previous study (cf. Figures 11-1 to 11-7 of Southworth and Sekanina, 1973), it 
is still too low to disclose the discrete levels that showed up in the photographic sam- 
ple. The resolution power has not been enhanced either by applying more restrictive 
criteria (e. g. , by accepting only well-observed, fully reduced meteors and/or those 
at small enough zenith distances), or by dividing the sample into physically more 
homogeneous groups (e.g. , by separating "bright" meteors from "faint" ones, etc.). 
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The only four features that we can detect on the height-velocity diagrams are those 
reported by us previously (Southworth and Sekanina, 1973): the separation between the 
major cluster of meteors at low V and the minor cluster at high V , the considerably 

OO OO 

larger scatter in the major cluster, the minor cluster more pronounced among the 
brighter meteors, and the absence of the high-velocity cluster for small -perihelion 
meteors. 


4. 3 Results from Least-Squares Solutions 


The negative results from the height-velocity plots led us to attack the problem 
numerically. We divided the whole sample into a number of groups and forced, for 
each group, a fit of the form 


h = hgQ + A log 



(4-1) 


where h_^ is the height at V =30 km sec In addition, we also made runs with the 

oU oo 2 

retained quadratic term [log (V /30)] and found that it was never important. 


OO 


The results of the linear solutions are included in Table 4-1. For each group, the 
table gives the number of meteors N, the height at V = 30 km sec ^ (^ 30 ^ mean 

error, the height slope A and its mean error, and the standard deviation of the individual 
meteor heights from the solution. 


Most of the data in the table refer to the diffusion height at maximum ionization, 
but five groups give solutions for the beginning height: Groups 2 , 4, 8 , 10 , and 12. 
These runs are marked by an asterisk in the table. The individual groups include: 

Groups 1 and 2. All meteors. 

Groups 3 and 4. Meteors for which geometric heights were found. Being more 
fully reduced, these meteors are considered to give better precision. 

Group 5. Meteors for which geometric heights were not found. The data for these 
meteors are considered, on an average, to be of lower precision. These meteors are 
not used in any of the groups after 6 . 
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Table 4-1. Results of the least-squares solutions of the height— velocity diagrams. 


Group 

number 

N 

^30 

A 

m. e. 

1 

14502 

89.87 ± 0.03 

19. 1 ± 0.2 

±3.6 

* 

2 

14502 

93.76 ± 0.03 

18.6 ±0.2 

±3.7 

3 

7639 

89.77 ± 0.04 

19.7 ± 0.3 

±3.3 

* 

4 

7639 

93.68 ± 0.04 

19.2 ± 0.3 

±3.5 

5 

6863 

90.01 ± 0.05 

18.4 ± 0.3 

±3.8 

6 

7299 

89.61 ± 0.04 

18.0 ± 0.3 

±3.6 

7 

3944 

89.59 ± 0.05 

19.0 ± 0.4 

±3.4 

* 

8 

3944 

93. 92 ± 0.06 

17.7 ± 0.4 

±3.5 

9 

4410 

90.21 ± 0.05 

19.5 ± 0.3 

±3.2 

10 

3226 

93. 59 ± 0.06 

19. 1 ± 0.4 

±3.4 

11 

2468 

90.08 ± 0.06 

19. 1 ± 0.4 

±3.3 

12 

1900 

93.66 ± 0.08 

18. 1 ± 0.5 

±3.5 

13 

2451 

90.46 ± 0.07 

15.6 ±0.5 

±3.2 

14 

1959 

90.38 ± 0.07 

24.3 ± 0.5 

±3. 1 

15 

746 

89.90 ± 0. 14 

20. 1 ± 1.3 

±3. 1 

16 

817 

90.31 ± 0. 11 

17. 9 ± 0.9 

±3. 1 

17 

844 

90.35 ± 0. 11 

20.0 ±0.8 

±3. 2 

18 

2003 

90.23 ± 0.07 

19. 9 ± 0.4 

±3.3 

19 

1876 

90.36 ± 0.07 

15.7 ± 0.5 

±3. 1 

20 

1030 

90.35 ± 0. 10 

21.0 ±0.7 

±3.2 

21 

1504 

90. 10 ± 0.08 

22.2 ±0.5 

±3.2 

22 

952 

90. 14 ± 0. 11 

22. 4 ± 0.6 

±3.2 

23 

596 

90. 15 ± 0. 14 

22. 5 ± 0.8 

±3.3 


* 

This run refers to the beginning height. 


Group 6. All meteors at zenith distances less than 45°. 

Groups 7 and 8. Well-reduced meteors at zenith distances less than 45°. 
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Groups 9 and 10. Well-reduced meteors with reliable ionization curves. 

Groups 11 and 12. Well-reduced meteors with reliable ionization cuiwes and at 
zenith distances less than 45°. 


The following groups are 
Group 13. Meteors with 
Group 14. Meteors with 
Group 15. Meteors with 
Group 16. Meteors with 
Group 17. Meteors with 
Group 18. Meteors with 
Group 19. Meteors with 
Group 20. Meteors with 
Group 21. Meteors with 
Group 22. Meteors with 
Group 23. Meteors with 


subdivisions of group 9: 

intrinsic magnitudes at maximum ionization M 11. 5. 
M > 11.5. 

perihelion distances q < 0. 25 a. u. 

0. 25 < q < 0. 50 a. u. 

0. 50 < q < 0. 75 a. u. 
q > 0. 75 a. u. 

aphelion distances q' s 2. 5 a. u. 

1. 5 < q' < 2. 5 a.u. 
q' < 1.5 a. u. 

q' < 1. 5 a.u. and q > 0. 5 a.u. 
q' < 1. 5 a. u. and q > 0. 75 a. u. 


We make the following conclusions from the results in Table 4-1: 

A. The beginning heights are about 3. 5 to 4 km above the heights at maximum 
ionization. (This difference is relatively smaller for meteors with reliable ionization 
curves only because of the reduction procedure, described in Section 6. 1. ) 

B. The meteors for which the geometric heights could be found appear to be of 
better quality, as expected. 

C. The meteors with reliable ionization curves appear to be only slightly better 
than those whose curves are less reliable. 

D. The zenith distance does not have any improving effect on the solution. We 
note, however, that it is difficult to reconcile numerically the effects of zenith-distance 
selection in Table 4-1 with the results in Tables 5-3 and 6-1; and we tentatively conclude 
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that the selection process for meteors in Table 4-1 was somehow anomalous with 
respect to zenith distance. 

E. A sharp difference occurs between the height slope A of the bright meteors 
and that of the faint ones. This effect is apparently due to the fact, already mentioned, 
that the high-velocity cluster of meteors is more pronounced among the brighter meteors. 
However, we have not yet taken the selection effects of diffusion and recombination into 
account in this part of the analysis. 

F. The height slope A is on an average 20 ± 5 and corresponds to the density- 

velocity dependence p ~ with n = 1. 5 ± 0. 4. This is considerably less than 

00 

Ceplecha's (1968) n = 2. 5. However, qualitative allowance for the selection effects of 
recombination and diffusion would result in closer agreement with Ceplecha. 

G. The meteors with q < 0. 25 a. u. appear to have somewhat smaller heights, but 
the difference from the other groups is only about 0. 4 km. 

H. The aphelion distance, like the perihelion distance, appears to have only a 
small effect on the heights. In both cases, however, the effects are formally signif- 
icant in a statistical sense, so that further study might well be repaid. 
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5. FRAGMENT ATION 


5. 1 Introduction 


The research described here constitutes a continuation of our previous studies of 
fragmentation, and certainly does not reach the end of the problem. The limitation of 
time unfortunately prohibited use of the measured amplitudes of individual Fresnel 
extrema, which are the best available data for fragmentation studies. Instead, we used 
the measured number of Fresnel extrema and have obtained the first set of values of the 
spread of meteoroid fragments along the trajectory, but mathematical peculiarities 
prevent our computing the complete distribution of values of the spread. 

5. 2 Inferring Fragment Spread from the Observed Number of Extrema 

This section describes computations of the number of Fresnel-pattem extrema 
that a human measurer or the pattern-measuring computer prc^ram would measure, 
taking into account diffusion and fragmentation in the ion trail and noise and limited 
d 3 mamic range in the radar equipment. A Fortran program was written for the com- 
putations, which were then to be used to infer the spread of the fragments along the 
meteor trajectory. 

McKinley (1961, p. 191) showed that when a meteor has passed the specular 
reflection point by a distance x, the signal amplitude (the power or the voltage squared) 
returned to the radar receiver is 

where 

F = (5-2) 
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is the effective length of the principal Fresnel zone, X is the radar wavelength, R is 

the range from the radar to the electron column, q is the number of electrons per unit 

length in the column, and C and S are Fresnel integrals. Here we have assumed that 

the electrons are uniformly distributed along the electron column within the interval 

2 

of interest (a few times F). The unit of is the response to a single electron at 
distance R. The subscript u denotes that diffusion and fragmentation have been 
neglected. Without loss of generality for this purpose, we shall assume 

q=l (5-3) 


since only relative values of A enter into determining the number of extrema. McKinley 
introduced the Cauchy approximations for the Fresnel integrals 


(^) 

Q - 1 ^ 

S — 1 cos 1 I 

Try V 2 / 


9 


9 


(5-4) 

(5-5) 


where 


y = ^ X . (5-6) 

2 

The successive maxima i = 1, 2, ... of A^ are found at approximately 

y^ = 4i-| . (5-7) 

Equations (5-4) and (5-5) are sufficiently accurate for the present purpose when y > 1, 
i. e. , at and after the first Fresnel maximum. 

Study of Loewenthal's (1956) diagrams shows that an alternative representation of 
A can be found in the vector sum of two components, ("deferent”) indicating the 
response of the principal Fresnel zone and ("epicycle") indicatii^ the fluctuating 
integral response of the later zones; E^ rotates continuously with respect to D^. 
Representing the relative phase by 4>, we have 
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(5-8) 




+ 2D E 
u u 


cos 4> . 


Substituting equations (5-2) throu^ (5-6) into (5-1) and comparing tenns with equation 
(5-8), we find 


D^=F 


F 

u 2irx 


and 



(5-9) 

(5-10) 


(5-11) 


Accurate integrations by Southworth (1962a, b), i. e. , without use of the Cauchy 
approximations, showed the validity of the epicycle— deferent representation; it will 
be equally valid to use the Cauchy approximations here because our applications are 
limited to observations of 5 to 30 extrema (3. 1 < y < 7.7), where the approximations 
are very good. 

Ambipolar difftisim of the electrons and ions away from the meteroid trajectory 
leads to attenuation of the radar return from the electron column. McKinley showed 
that the return at time t after formation of the column is reduced by a factor 


a (t) = exp 


32 TT^ Gt\ 


(5-12) 


where G is the ambipolar diffusion coefficient. Greenhow and Neufeld's (1955) obser- 
vational representation of G is 

log^Q G (cm^ sec"^) = 0. 068 h - 1. 67 , (5-13) 


where h is the height above sea level (km). Equation (5-13) is reasonably good for 
synoptic-year meteors and is exact by definition for the 1961 — 65 meteors, where h was 
determined from G and equation (5-13), 
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We take diffusion into account in this analysis by introducing appropriate attenua- 
tion factors in the D- and E -components of the returned signal. The D-component is 
affected by diffusion during the time t = x/V (where V is the meteoroid velocity) since 
the meteoroid passed the specular reflection point; accordingly, we can multiply D by 
the factor 


a = exp 


16 TT-^ Gx 1 


(5-14) 


Introducing Loewenthal's (1956) parameter 


C = 


SttGF 


(5-15) 


we then have 


a^=exp(^-^^) . (5-16) 

The E -component is affected by diffusion over a much shorter time, and accurate 
calculation of the relatively small effect is complicated. For completeness, however, 
we estimate that the E-component is effectively formed one-quarter revolution earlier; 
by differentiation of equation (5-11), this corresponds to a distance 



earlier; and by comparison with equation (5-16), we multiply E by the factor 


c = exp 


ttCF 

x 


(5-18) 


Fragmentation of the meteoroid and the consequent spread of the fragments along 
the trajectory by differential deceleration replace the Fresnel pattern of a single 
meteoroid with the sum of the Fresnel patterns of the separate fragments relatively 
out of phase. The effect is to reduce the E -component of the returned signal. We 
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model the fragment spread by a Gaussian distribution of fragments about the fragment 
mean, with a root-mean-square (rms) deviation a from the mean. Then, the relative 
number of fragments s ± ds from the mean is 



The relative phase of a fragment at x + s, with respect to the fragment mean at x, is 


[(x+s)2-x^] 

[cf. eq. (5-11)], which we approximate by 
, _ 2ir sx 

Then the relative amplitude of the E -component is given by 

CO 


I 


n(s) cos <|)g^ ds 


b = 


oo 

/ 


= exp 


n(s) ds 



“OO 


(5-20) 


(5-21) 


(5-22) 


and we can accordingly multiply E by b. 

In theoretical Fresnel patterns (with the use of the above assumptions), the 
oscillating component never vanishes in the direction of increasing x, but there is 
nonetheless a definite bound to the number of observable extrema. Beyond that bound, 
the slope of the returned signal, though oscillating, is always negative. Both the 
human film measurers and the computer program that later measured Fresnel patterns 
started at the first maximum and proceeded as far as possible. Both stopped when 
there were no further extrema, as just described, even though an oscillating slope 
could often still be seen. They also stopped when there were too few pulses per 
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extremum or too much noise in the data to recognize the next extremum, as well as 
when the returned signal became too small to be recorded. Moreover, the human 
measurers (1961-65 meteors) stopped at 20 extrema, and the program (synoptic-year 
meteors), at 30. 

To find the theoretical number of extrema, we computed the slope of the Fresnel 
pattern at points where c() = 3 tt/2 (modulo 2 -rr), i.e. , approximately halfway between 
each minimum and the following maximum. This is the point where the slope has a 
local maximum in the Fresnel pattern of a meteor unaffected by diffusion and fragmen- 
tation and approximately where it has a local maximum when diffusion and fragmenta- 
tion are brought into account. Denoting these points by m = 1, 2, . . . , according to 

th 

the order in which they follow the m''^^ minimum, we have 


and 


(j) = TT (2m - 
x=F(2m+^) 


(5-23) 


(5-24) 


Introduction of a [eq. (5-14)], c [eq. (5-18)], and b [eq. (5-22)] transforms equation 
(5-8) to 


= (aD^)^ + (bcE^)^ + 2 abcD^E^ cos cj) . (5-25) 


Substitution of equations (5-9), (5-10), and (5-11), differentiation, and substitution of 
equation (5-24) yield the approximate local maximum slope 



2F 


da (c/27t)^ 
dx 2 m + 1/4 


db abc 
dx F 


(bc/2TT)^ 

F(2m+1/4)^^^J 


(5-26) 


where dc/clx has been neglected. 

The accuracy with which the above slope could be measured was estimated as 

2 2 

follows. We are, in fact, concerned with the relative slope (1/A ) [d(A )/dx], i.e.. 
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2 

the logarithmic derivative of A , because all signals were recorded in logarithmic 
units and the "noise" of recorded signals was more nearly constant in a logarithmic 
scale than in an absolute scale. The slope in question occurs near a point of inflection 
of the Fresnel pattern, where the observed pattern is sensibly straight over an inter- 
val of roughly one-sixth the local period Ax of the Fresnel oscillations. To the 
human measurer, the uncertainty e in measuring this slope is thus approximately 
the xmcertainty e^ in the difference of the signal at the two ends of this interval, 
divided by the length of the interval, or 


6e 


d 


e = A 
s Ax 


^®d 

F [(2m + 5/4) - (2m - 3/4) 


(5-27) 


where we have used equation (5-24) to evaluate the length of the interval. We estimate 
e^ to be 


®d = 


0. 05 N 


(P-3) 


V2 


(5-28) 


2 

where 0. 05 is the relative change in A corresponding to a single step of the digital 
recording equipment near the middle of the dynamic range of the receivers, N is the 
noise of a single recorded signal in digital-recording-equipment units, and p is the 
number of pulses between successive maxima. Our estimate of the effective value of 
N for the synoptic year is 1, and the equivalent value for 1961 — 65 (when digital 
recordii^ was not used) is 2. Nonetheless, N is left as a free parameter in the com- 
putations. Since 738 pulses sec”"^ were transmitted, we have 


= 738 F [(2m +5/4)^^^ - (2m -3/4)^^^] 
P V 


(5-29) 


Both measuring techniques considered approximately one cycle of the Fresnel oscilla- 
tion in measuring each extremum, and both were instnacted to stop measuring if there 
were fewer than 3 pulses per cycle. Thus, p - 3 is appropriate in equation (5-28); the 
1/2 power is used in the expectation that the mean error is inversely proportional to 
the square root of the number of observations. 
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The computer program to determine the observable number of extrema evaluated 


for successive values of m = 1, 2, . . . . A value T(m) = 0 implies that 

n = 2m + 1 (5-31) 

extrema could be observed, to wit the m maxima and m minima preceding the point 
at which T(m) = 0 and the following maximum. Strictly speaking, a value T(m) > 0 
followed by T(m + 1) < 0 also implies that n = 2m + 1 extrema are observable. How- 
ever, the subsequent use of n requires that it be a continuous variable, so that n is 
determined in that case by inverse linear interpolation as 


n = 2m + 1 + 


2 T(m) 

T(m) - T(m+ 1) 


(5-32) 


An approximate correction for the limited d 3 mamic range of the receivers was 
introduced. We are concerned only with the lower end of the range, because echoes 
bounded at the upper end were not reduced. The number of meteors increases rapidly 
with decreasing mass, so that most of the meteor echoes are near the faint observ- 
ing limit. We estimated that the average reduced meteor was L = 10 db above this 
limit at its first maximum. It was not possible to test the estimate, as it would entail 
an elaborate analysis of calibration data, but trial computations showed that computed 
values of n were not at all sensitive to the estimate. The correction consisted of 
subtracting from T(m), when A /F < 10 ' , the quantity 

T' = 0. 1 log^Q (A“^ 10^^^) , (5-33) 

which simulates the added difficulty of observing extrema after an echo has decayed 
by more than 10 db. 
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The procedure described above determines the number n of observable extrema 
as a function of five parameters: C, ff/F, F/V, N, and L; the computer program 
created a three-dimensional table of n as a function of C, a/F, and F/V for given 
values of N and L. The left-hand portion of Table 5-1, labeled "Extrema, " is a 
severely abbreviated example of the results. The parameter C is dimensionless; C = 0 
denotes no diffusion, and C > 0. 3 denotes so much diffusion that successful reduction is 
not likely. The left-hand column, "Frag.," designates cr/F, where F ranges from 0. 5 
to 1. 2 km; thus we are concerned here with rms fragment distances, from the fragment 
mean, of up to a few tenths of a kilometer. Table 5-la, F/V = 0. 015, is relevant 
to fast meteors (velocity V = 33 to 80 km sec~^), and the other sections, to slower 
meteors. Values of n > 100 have been truncated to 100, and those for n < 3 were 
computed by using a conventional value for T(0) = 0.48. The letter R denotes cases 
where the E -component exceeds the D-component; this is important in determining 
winds, but it is not important for the present purpose. 

Examination of Table 5-1 shows that n is not always a monotonically decreasing 
fvinction of a/F, but sometimes is of the form shown in Figure 5-1. When, as at Q 
in Figure 5-1, n has a relative minimum n^, a/F can be determined from n only if 
n < n^. When n > n^, all that can be said is that a/F < s^. The right-hand half of 
Table 5-1, labeled "Fragment Spread, " gives values of a/F as a function of n, obtained 
by inverse interpolation of an imabbreviated version of the entries under "Extrema. " 
When n has a relative minimum, as in Figure 5-1, the table uses P Q R S for computing 
a/F, but only the section R S will later be used. 

5. 3 Results on Fragment Spread 

We used the procedure just described to compute, as far as possible, a value of 
the fragment spread a for each of 17, 917 synoptic-year meteors and for each of 9691 
meteors observed in 1961 — 65 that had meaningful diffusion heights. The synoptic-year 
sample is, of course, less biased by observational limitations, but the 1961 — 65 sample 
represents larger meteors and is also more homogeneous in heights. 


47 


Table 5-1. Computed numbers of extrema. 
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Figure 5-1. An example of the computed observable number n of extrema as a function 
of fragment spread a in units of the effective length F of the principal 
Fresnel zone. See text for discussion. 

For each meteor, we used the mean number of extrema measured among all 
Fresnel patterns accepted by the reduction programs; thus we obtain only a mean meas- 
ure of the fragment spread. Unfortunately, it was not possible to use the numbers of 
extrema at individual stations and, therefore, not possible to determine changes in the 
spread. Moreover, no Fresnel pattern with less than five extrema was ever reduced, 
so there is a bias against our finding lai^e values of a. 

In view of the significant proportion of meteors for which u could not be computed, 
we carried out analyses in two forms. First, where c was determined, we attempted 
to find its functional dependence, if any, on other variables. Second, dividing the 
meteors into three classes — 1) those with a greater than some given value, 2) those 
with O’ smaller than the given value, and 3) those where o could not be determined — 
we attempted to see if other properties showed any meaningful difference among the 
classes. All in all, a very extensive set of statistics was compiled, but almost all 
of them revealed nothing, at least to us. In particular, no important relationships 
seem to exist between o and perihelion distance, aphelion distance, or inclination. 
However, the data should be reexamined when we understand better the results that 
have been found. 
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Table 5-2 contains the coefficients found by least-squares fits of log^Q a to func- 
tions of velocity, radar magnitude, radiant zenith distance, and mean height above 
sea level. We do not interpret them, and they constitute data for a future theory of 
fragmentation. Values of cr < 0. 01 were omitted from the analysis because they are 
much too sensitive to observational errors. Undetermined and negative values were, 
of course, also omitted. Other values of noise give qualitatively similar results. 

The relatively larger errors of a single observation ("s. e. 1") for the synoptic year 
probably reflect inhomogeneities in the height data used, leading to errors in the value 
of C inferred from the heights and thence in the values of a inferred from C. 

Table 5-2. Least-squares fits for cr (km). 


Coefficient of 


Data 

Number of 
meteors 

Meteors 

usable 

1 

log^(j(V/30) 

M - 11 logjQ(cos Z/0.707) h - 90 

s. e. 1 

1961-65 

9690 

8062 

-1. 08 

-0.20 

-0. 01 

+0.05 

+0. 007 

0. 17 

noise = 2 



±0.03 

±0. 02 

±0.00 

±0.01 

±0. 000 


1961-65 

9690 

6260 

-1.20 

-0. 36 

-0.02 

+0.02 

+0.010 

0.21 

noise = 4 



±0. 05 

±0. 02 

±0. 00 

±0.02 

±0. 001 


Synoptic year 

17576 

10909 

-0. 94 

-0.31 

-0. 02 

-0.00 

+0. 005 

0.30 

noise = 1 



±0. 04 

±0. 02 

±0. 00 

±0.02 

±0. 000 


Synoptic year 

17576 

10017 

-1. 05 

-0. 37 

-0.02 

-0.00 

+0.010 

0.29 

noise = 2 



±0.02 

±0.02 

±0. 00 

±0.02 

±0. 000 



Table 5-3 contains coefficients of least-squares fits of mean height to the other 
variables used in Table 5-2. It was prepared for another purpose, but is presented 
here as further data for a future theory. 


Table 5-3. 


Least-squares fits for-^ (h^ 


eg 


+ h 


end 


) (km). 


Coefficient of 


Data 

Number of 
meteors 

1 

logio(V/30) 

M - 11 

logio(cos Z/0.707) 

s. e. 1 

1961-65 

9691 

+ 87. 04 

+22. 18 

-1. 94 

-3.43 

6. 02 



±0. 19 

±0.46 

±0. 08 

±0.47 


Synoptic year 

17917 

+ 90. 84 

+19. 96 

-0.38 

-5.08 

7. 16 



±0. 06 

±0.37 

±0. 06 

±0.41 
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Table 5-4 presents an example of the second form of the statistics described 
above; moreover, it contains a most puzzling result. On any simple model, we should 
expect a smooth inverse relation between velocity and spread, but here we observe 
that spread is largest for 16 V < 25, next largest for 40 ^ V, next for 25 ^ V < 40, 
and smallest for V < 16. This vmexpected relationship provides not only more data for 
a future theory, but also an obstacle to further analysis until it is elucidated. Thus 
far, we have drawn only two minor conclusions from Table 5-4: first, that the con- 
fusing appearance of many of the statistics we computed and then did present here 
may be explained by the unexpected relation between V and a; second, that our previous 
estimates of N = 2 for 1961 — 65 and N = 1 for the synoptic year are better than N = 4 
and N = 2, respectively, as evidenced by the large proportion of undetermined a's at 
the large values of N. 

Table 5-4. Relative numbers of meteors having a ^ 0. 1 km. 






1961- 

65 


Synoptic year 




Noise = 2 


Noise = 4 

Noise = 

1 Noise = 2 



a < 0. 1 

0.81 


0. 84 

0. 85 

0. 90 


V < 16 

CT > 0. 1 

0. 19 


0. 11 

0. 15 

0. 09 



undetermined 

0.00 


0.05 

0. 00 

0.01 



n 


151 



951 



a < 0. 1 

0. 56 


0.58 

0.40 

0.49 

16 < 

V < 25 

a > 0. 1 
undetermined 

0.43 


0.31 

0. 58 

0.47 

0.01 


0. 11 

0.02 

0.04 



n 


1506 



4467 



a < 0. 1 

0. 64 


0. 53 

0.77 

0. 60 

25 £ 

V < 40 

a > 0. 1 
undetermined 

0.28 


0. 11 

0. 19 

0. 12 

0.08 


0.36 

0.04 

0. 28 



n 


5817 



8963 



a < 0. 1 

0. 48 


0.45 

0. 58 

0.61 

VI 

o 

V 

a > 0. 1 

0.36 


0. 14 

0,33 

0. 12 



undetermined 

0. 16 


0.41 

0.09 

0. 27 



n 


2216 



3195 
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Table 5-5 contains data comparable to Table 5-4, except that the dividing value 
for a is 0. 03 km. It does not contain a surprise as Table 5-4 did, but shows that very 
small values of a are relatively xmcommon. 


Table 5-5. Relative numbers of meteors having a ^ 0. 03 km. 




1961-65 

S3moptic year 



noise = 2 

noise = 1 


CT < 0. 03 

0.00 

0.00 

V < 16 

a > 0. 03 

1.00 

1.00 


undetermined 

0. 00 

0.00 


n 

151 

951 


a < 0. 03 

0.02 

0. 02 

16 £ V < 25 

cr > 0. 03 

0. 97 

0.94 


undetermined 

0.01 

0.04 


n 

1506 

4467 


(7 < 0. 03 

0. 05 

0. 03 

25 < V < 40 

a > 0. 03 

0.82 

0.38 


undetermined 

0. 13 

0. 59 


n 

5817 

8963 


CT < 0. 03 

0.07 

0. 05 

40 < V 

CT > 0. 03 

0.83 

0. 78 


undetermined 

0.10 

0. 17 


n 

2216 

3195 


5.4 Conclusions 


We badly need a physical theory of fragmentation to interpret the observations 
presented in Section 5. 3 and thus to understand other observations. The amplitudes 
of the Fresnel patterns from the synoptic year should be analyzed. The approximate 
magnitude of the fragment spread along the trajectory found in our previous report 
(Southworth and Sekanina, 1974, Table 2) is well confirmed. [Note the differing defini- 
tions of w (1974) and a (this report). ] 
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The conclusion in our 1974 report that we have no evidence for observational 
selection by fragmentation is shakily confirmed. That conclusion rested on the bunch- 
ing of the data to low values of the spread and, in particular, on the absence of 
higher spreads at low velocity. The analysis reported here shows a broader distribu- 
tion of the spread than was envisaged before. Moreover, while few large spreads 
occur at the lowest velocities, many do in the interval 16 < V < 25. 
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6. OBSEEVATIONAL BIASES 


6. 1 Radiant Zenith Distance 


Southworth and Sekanina (1974, p. 12) included a note added in proof stating that 
we had just realized a systematic bias due to the fact that the radar system was not 
oriented to observe long trails from meteors with small cos Z^^. The orientation was, 
of course, in the original design of the system, to enable it to observe all of some 
long trails with a minimum n\imber of stations. 

In order to determine the necessary correction, we computed, for each meteor 
in the synoptic year, the projected extent of the station array on the meteor trajectory. 
When a perpendicular is dropped onto the trajectory from each station, the projected 
extent is the distance between the most widely separated perpendiculars. The specular 
reflection point for any station is approximately halfway between the perpendicular 
from that station and the perpendicular from the transmitter station; accordii^ly, the 
maximum distance between specular reflection points is approximately half the pro- 
jected extent of the station array. To find the ionization curve, the computer program 
fitted a quadratic in distance along the trajectory to the logarithm of the line density 
determined at each specular reflection point. If this quadratic was concave downward, 
and if it intersected the minimum detectable line density no more than 4 km before the 
first specular reflection point and again no more than 2 km after the last observed 
Fresnel extremum, the ionization cuiwe was classified as well-determined. A quad- 
ratic intersecting the m inimum detectable line density outside either limit was trun- 
cated at that limit, and the curve was classified as not well-determined. A quadratic 
concave upward was replaced by a rectangle bounded at those limits and classified as 
not well-determined. 

The trail length is the distance between the intersections of the quadratic with the 
minimum detectable line density or the ends of the curve as limited above. Allowing 
roughly 2 km for the distance between the last specular reflection point and the last 
observed Fresnel extremum, a not -well-determined trail length is therefore limited 
to, at most. 
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Li ^ 1 , “ 0. 5 E„ + 8 km , 
length p ’ 


(6-1) 


where E is the projected extent of the station array. Well-determined lengths are 
P 

systematically shorter than those that are not. 

Some exceptions arise to the limit expressed in equation (6-1). In an attempt to 
correct for minor observational deficiencies, those ionization curves that had been 
truncated at either the beginning or the end (not both) by less than half the distance 
between the maximum and the intersection with the minimum detectable line density 
had that truncation restored for this purpose, and the curve was reclassified as well- 
determined. 


We divided the meteors into groups; E^ < 10 km, 10 < Ep < 20, 20 < Ep < 30, 
and 30 < Ep; we also divided them by magnitudes ^ 11 and by well -determined and not- 
well-determined ionization curves. In each group, we then fitted coefficients A and 
B by least-squares of the form 


fcos Z 


R 


■^^\eg~\nd^ ^ ^ ® ^’^\0.7071 


(6-2) 


where hj^^^ and h^^^ are the heights above sea level of the beginning and end of the 
trail. Table 6-1 gives the values foimd for A and B and various other data. Table 6-la 
gives well-determined ionization curves, and Table 6-lb presents those that are not 
well-determined. The logarithmic (geometric) mean is denoted by L MN, and the anti- 
logarithm of the rms deviation of the logarithm from the mean of the logarithms, by 
L SD. MAX and MEST are the largest and smallest values in the group; PR. 1 is the 
error that will be exceeded in 10% of trails, using Student's criterion; and PROJ is 
the mean projected station extent. The numbers of meteors with ionization curves 
extended for the purpose of this study, and therefore not strictly subject to the limit 
in equation (6-1), are designated (BEG) and (END). 


Figure 6-1 shows the regression lines for well-determined ionization curves 
fainter than magnitude 11. They have been terminated at the 1/e points in the distri- 
bution of in (cos Z /0. 7071) to illustrate the range of cos Z„ occurring In each group. 

XV 
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Table 6-1. Statistics of vertical trail length and cos 


a) Well-determined ionization curves 



GOOD ION C 

GOOD ION C 

GOOD ION C 

GOOD ION C 


PROJ 

GT 30 

PROJ 20-30 

PROJ 10-20 

PROJ 

LE 10 

MAGN 

GT 11 

LE 11 

GT 11 

LE 11 

GT 11 

LE 11 

GT 11 

LE 11 

NOBS 

273 

163 

3641 

2036 

1382 

851 

86 

60 

(BEG) 

43 

16 

574 

206 

204 

91 

4 

0 

(END) 

10 

10 

248 

176 

103 

70 

6 

0 

HB-HE 

L MN 

7.81 

8.65 

8.29 

8.93 

6.43 

6. 98 

7. 17 

6.66 

L SD 

3.00 

3.21 

3. 10 

3.25 

2.70 

2. 89 

1.70 

1. 19 

MAX 

19.4 

18.8 

23.0 

25.6 

19.6 

20.2 

14.2 

10.0 

MIN 

3.0 

3. 1 

1.2 

2.2 

0.9 

0. 9 

4. 1 

4.7 

COSZR 

L MN 

0.655 

0.670 

0.715 

0.721 

0. 587 

0.632 

0.796 

0.818 

L SD 

0. 103 

0. 091 

0. 142 

0. 135 

0. 230 

0.238 

0. 101 

0.091 

MAX 

0.75 

0.75 

0.90 

0.90 

0. 94 

0. 94 

0. 93 

0, 93 

MIN 

0.29 

0.31 

0. 12 

0. 18 

0.09 

0. 10 

0. 54 

0.58 

A 

2. 1120 

2.2002 

2. 1072 

2. 1799 

2. 0044 

2. 0298 

1. 8909 

1. 8694 

PR. 1 

0. 0407 

0. 0496 

0. 0095 

0.0129 

0.0143 

0.0173 

0. 0546 

0. 0639 

B 

0. 7408 

0.7889 

0.6919 

0.4748 

0.7721 

0. 7725 

0. 6730 

0. 1874 

PR. 1 

0. 2343 

0. 3407 

0. 0478 

0. 0686 

0. 0330 

0. 0440 

0.3154 

0. 3507 

PROJ 

31.45 

31.37 

25. 06 

24.91 

15.78 

15.80 

8. 17 

8.21 


LOWBND I C 

b) Not well-determined ionization curves 
LOWBND I C LOWBND I C 

LOWBND I C 


PROJ GT 30 

PROJ 20 

-30 

PROJ 10 

-20 

PROJ LE 10 

MAGN 

GT 11 

LE 11 

GT 11 

LE 11 

GT 11 

LE 11 

GT 11 

LE 11 

NOBS 

352 

97 

3674 

1361 

2701 

983 

134 

68 

(BEG) 

0 

0 

0 

0 

0 

0 

0 

0 

(END) 

0 

0 

0 

0 

0 

0 

0 

0 

HB-HE 

L MN 

8.36 

9.21 

9.07 

9. 87 

6.75 

7.26 

8. 51 

9. 06 

L SD 

2.46 

3.31 

2. 55 

2.71 

2. 55 

2. 90 

1.30 

1.87 

MAX 

19.2 

19. 1 

20.9 

19.2 

23.0 

20.6 

16.7 

17.8 

MIN 

1.7 

3.5 

1.2 

1.7 

1.2 

1.4 

6. 1 

5.8 

COSZR 

LMN 

0.656 

0.637 

0.698 

0.716 

0. 548 

0. 563 

0.794 

0.812 

L SD 

0.112 

0. 117 

0. 155 

0. 141 

0.215 

0.227 

0. 100 

0. 100 

MAX 

0.75 

0.75 

0.90 

0.90 

0.94 

0.94 

0.94 

0.94 

MIN 

0.20 

0.33 

0. 12 

0. 11 

0.09 

0. 10 

0. 53 

0. 53 

A 

2. 1897 

2.3056 

2.2149 

2.2808 

2. 1259 

2. 1820 

2.0875 

2.0866 

PR. 1 

0. 0242 

0. 0637 

0.0061 

0.0105 

0. 0067 

0.0112 

0. 0278 

0. 0544 

B 

0. 8899 

0.8195 

0.7635 

0.7149 

0. 8507 

0. 8758 

0.4610 

0. 8499 

PR. 1 

0. 1300 

0.3035 

0. 0275 

0. 0536 

0.0142 

0. 0242 

0. 1627 

0.2952 

PROJ 

31. 53 

31.55 

24.65 

24.51 

15.41 

15. 74 

8. 57 

8.03 
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Lines for brighter meteors, or for not -well-determined ionization curves, are similar, 
but they are systematically displaced upward and are rather less regular. It is evident 
that a regression line fitted to all the data would show a steeper slope and that it would 
not properly represent the physical relationship. Moreover, there is still reason to 
fear that some of the separate regression lines are also subject to selection effects. 
Dotted lines in Figure 6-1 show the upper bound within the figure for a meteor 

of each group that has the average projected station extent for that group and that did 
not have its ionization curve arbitrarily extended, which is 


It seems likely that the regression line for the group 10 < E < 20 has been biased. 

In order to correct our statistics, we must decide on a relation between trail 
length and cos Z_. Taking an approximate weighted average of the values of B in 

JA 

Table 6-1 with E^ > 20, we adopt 

B = 0. 7 ± 0. 1 . (6-4) 


6. 2 Projected Station Extent 

We need further to correct for the tendency of the limitation in projected station 
extent to cause us to underestimate the trail length and thus the mass of a meteor. 

(It is not, however, necessary to Include any new correction for the size or shape of 
the atmospheric volume obseiwed by the radar system. ) Figure 6-2 is a plot of the 
A-coefficients in Table 6-1, and empirical straight lines have been drawn through 
points for well-determined and not-well-determined curves. In principle, both curves 
shovild have horizontal asymptotes at the right, but we accept them as practical 
approximations. The convergence of the lines definitely suggests that not-well-deter- 
mined curves would be well-determined if only the station array were large enough. 

To correct our statistics, we adopt the average slope 0. 18 of the two lines and esti- 
mate that the trail length is independent of E above E = 40. The corresponding 
correction factor for the system sensitivity that must be inserted into existing com- 
puter programs is 
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Figure 6-2. Coefficients A [see eq. (6-2) and Table 6-1] as functions of mean projected station extent E 



and the revised mean vertical trail length is 


\eg->'end=‘2.S(cosz/-'' . (6-6) 

Equation (6-6) replaces equation (16) from Southworth and Sekanina (1974), which we 
recc^nized as suspect while reading the proofs of the report, and which read 

V‘end='«-’(‘=<>= V"'®® • 

The differences, though real, are much too small to invalidate any of the conclusions 
we drew from it. 

6. 3 Antenna Sensitivity 

In Southworth and Sekanina (1973), we described a revision of the previously 
adopted radiant sensitivity pattern, based on measurements of antenna sensitivity. 
There is, however, a companion factor to the radiant sensitivity, which takes into 
account, as a function of declination, the interval of right ascension within which the 
radiant sensitivity is large enough to be applied, since radiants outside that interval 
are assigned zero weight. We find that this companion factor was not changed at that 
time; consequently, we have now changed it. The effect is not large. 



7. ORBITAL DISTRIBUTIONS 


7. 1 Data Samples 

This section presents tables of the distributions of velocity, radiant, and orbital 
elements of meteors observed in the synoptic year and in 1961-65. Most of the tables 
represent the data in three forms: 1) unweighted, 2) reduced to equal masses in the 
atmosphere, and 3) reduced to equal masses in space at the earth’s distance from the 
sun. We use three forms to facilitate comparison with the great variety of other 
published forms. Observations from earth satellites should be compared with the 
second form, and observations not restricted to the ecliptic (e. g. , comets), with the 
third form. 

The weights are computed as described in Southworth and Sekanina (1973), except 
for the changes discussed in Section 6. There is no correction for fragmentation. 

The correction for diffusion, empirically developed for the synoptic-year meteors, 
can be only approximate for the 1961—65 meteors. No correction could be intelligently 
made for the known bias against low velocities in 1961 — 65 caused by recombination 
(Southworth, 1973). 

Table 7-1 contains the number of observations used for the different samples 
presented. Meteors with slant ranges over 400 km were omitted from all samples 
because they are especially subject to errors caused by wind shears. Those with 
slant ranges under 70 km were not included, because such ranges are likely to be 
recording errors. Meteors outside the main lobe of the antenna pattern were left out 
because they have less accurate reductions and are likely to be unrepresentatively 
larger than meteors in the main lobe. A few 1961-65 meteors with seemingly 
erroneous masses were also omitted. Finally, we did not include any meteors from the 
space samples with radiants south of the ecliptic, because we cannot get a complete 
sample of southern radiants. 
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Table 7-1. Number of meteors used. 



S 3 moptic year 

1961-65 

Total reduced 

19698 

19180 

Observed sample 

14076 

14520 

Atmospheric sample 

14076 

14520 

Space sample 

12321 

11920 


Table 7-2 contains the unweighted distribution of mass in the S 5 moptic-year and 
1961-65 observed samples. The mean difference in equipment sensitivity between 
the two periods was approximately 22 db, corresponding to a factor of 12. 6 in mass. 
This factor, however, is far from being completely reflected in the mass distributions, 
because of the relative lack of low-velocity meteors in the earlier sample. 


Table 7-2. Distribution of mass. 


Mass 

(g) 

Synoptic year 

1961-65 




10 

0. 006 

0.048 

10~^ 

0. 096 

0. 326 

10"^ 

0.384 

0.457 

10~^ 

0.397 

0. 152 

10“^ 

0. 109 

0.019 

10"® 

0. 009 

0.000 

10"^ 




For the synoptic year, the tables in this section supersede the corresponding 
tables and figures in Southworth and Sekanina (1973) (the 1961—65 meteors have not 
been previously published in this form). Moreover, the space densities in the 1973 
report also require revision, which we hope to publish separately. 

7. 2 Velocities 


Table 7-3 contains the distribution of no-atmosphere velocities for the six samples. 
Each entry is the fraction of one sample having velocities in a 5-km sec ^ interval. 
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except that the last interval includes all velocities over 70 km sec Entries smaller 

than 0. 010 are in a compressed form of scientific notation; for example, 39-4 means 
-4 

39. X 10 (note the position of the decimal point). 


Table 7-3. Distribution of no-atmosphere velocity. 


Velocity 
(km sec”"^) 


l^moptic year 



1961-65 


Observed 

Atmospheric 

Space 

Observed 

Atmospheric 

Space 

10 

0. 029 

0.701 

0.256 

81-4 

0. 593 

0. 150 

15 

0.084 

0.214 

0.398 

0.037 

0.232 

0.323 

20 

0. 160 

0. 062 

0.221 

0. 101 

0. Ill 

0.281 

25 

0. 195 

0.017 

0.083 

0. 182 

0.043 

0. 148 

30 

0.201 

48-4 

0. 031 

0.225 

0.015 

0.068 

35 

0. 144 

12-4 

85-4 

0. 185 

42-4 

0.024 

40 

0. 068 

22-5 

18-4 

0. 099 

92-5 

54-4 

45 

0. 031 

47-6 

32-5 

0.046 

19-5 

11-4 

50 

0. 033 

23-6 

13-5 

0.040 

79-6 

25-5 

55 

0. 031 

13-6 

50-6 

0. 042 

49-6 

12-5 

60 

0.017 

43-7 

12-6 

0. 026 

20-6 

42-6 

65 

50-4 

87-8 

20-7 

79-4 

40-7 

38-7 

70 

50-5 

79-9 

14-8 

12-4 

35-8 

36-8 


The synoptic-year distribution is little changed from our 1973 report, with a then- 
unprecedented preponderance of low velocities; 1961 — 65 is closer to other observers* 
distribution. 


7. 3 Radiants 


Tables 7-4 through 7-9 show the radiant distribution in the six samples as a function 
of celestial longitude measured from the sun LAMBDA and celestial latitude BETA. 

Each entry is the mean radiant density, per unit area, in a 5® interval of latitude and 
longitude. The xmit of density is one-tenth the mean radiant density in the north 
celestial hemisphere; thus, "10" represents mean density. The entries have been 
smoothed with a triangular weighting function in longitude and latitude. The latitudinal 
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Table 7-4. Radiant distribution from the synoptic-year observed sample. 
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Table 7-4 (Cont. ) 
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Table 7-5. Radiant distribution from the 1961 — 65 observed sample. 
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Table 7-5 (Cent. ) 
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Table 7-6. Radiant distribution from the s 5 moptic-year atmospheric sample. 
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Table 7-6 (Cent. ) 
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Table 7-7. Radiant distribution from the 1961—65 atmospheric sample 
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Table 7-7 (Cent. ) 
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Table 7-8. Radiant distribution from the synoptic-year space sample. 
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Table 7-9. Radiant distribution from the 1961 — 65 space sample, 
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smoothing was continuous over the north celestial pole, and the longitudinal smoothing 
involved nearly equal great-circle intervals of longitude at all latitudes. In smoothing 
the "space" samples, which are truncated at the ecliptic, we assumed a mirror-image 
distribution in the southern hemisphere. 

Comparison of these tables with our 1973 report shows no large changes. (Note 
that the table entries in 1973 were differently defined, and harder to interpret. ) The 
1961-65 data are remarkably similar to those from the synoptic year. The strong 
concentration of relative velocities into two clusters (four, including the southern 
hemisphere) toward and away from the sun in longitude, but at high angles to the ecliptic, 
needs further discussion. 

7. 4 Orbital Elements 


Tables 7-10 through 7-15 present distributions of 1/a, e, q, q', i, and w, 
respectively, for the six samples. The first five show little difference from our 1973 
report for the synoptic year. Differences between 1961 — 65 and the synoptic year are 
evidently associated with the relative lack of slow meteors in the former period, but 
they have not been otherwise analyzed. 

Table 7-15 is new and is introduced to help imderstand Tables 7-6 and 7-9 (the 
space-sample radiant distributions). The space samples in Table 7-15 show a marked 
dip in the frequency of co near 0° and 180° ; in other words, it shows a tendency for the 
nodes to avoid the neighborhood of perihelion and aphelion. This is equivalent to 
avoiding \ - \q = 90° or 270° for the radiants. Qualitative explanations can be framed 
in terms of either Jupiter perturbations of parent comets or collisions near the 
ecliptic plane (for example, the asteroid belt), but they have not been worked out in any 
detail. 

Tables 7-16 through 7-21 give two-dimensional distributions of 1/a and e, 1/a and i, 
and e and i, for both space samples. They are independently scaled to a largest entry 
of 1000, but the row and column sirnis make it easy to change scales. In Tables 7-16 
and 7-17, the areas above the triangular array correspond to orbits with perihelion 
outside the earth's orbit, and the areas below, to orbits with aphelion inside the earth's 
orbit. The apparently hyperbolic orbits (e > 1) were probably displaced from the top 
right comer of the triangular array by observational error. 
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Table 7-10. Distribution of inverse semimajor axis 


1/a 


Synoptic year 



1961-65 


Observed 

Atmospheric 

Space 

Observed 

Atmospheric 

Space 

0.0 

0.010 

52-5 

13-4 

0.021 

94-5 

28-4 

0.1 

0.015 

25-4 

24-4 

0.021 

13-4 

31-4 

0.2 

0.028 

31-4 

71-4 

0.043 

77-4 

0.015 

0.3 

0. 056 

0.011 

0. 022 

0. 064 

0.015 

0. 033 

0.4 

0.085 

0.042 

0.060 

0.087 

0.027 

0.049 

0.5 

0. 102 

0. 051 

0. 083 

0.102 

0.073 

0. 104 

0.6 

0. 092 

0. 060 

0.088 

0.092 

0.098 

0. 107 

0.7 

0. 085 

0. 102 

0. 098 

0.084 

0.076 

0. 101 

0.8 

0.078 

0. 122 

0.113 

0.081 

0. 136 

0. 112 

0.9 

0. 080 

0. 153 

0. 127 

0. 079 

0. 171 

0. 127 

1.0 

0. 088 

0. 187 

0.133 

0. 080 

0. 174 

0. 124 

1.1 

0.088 

0. 140 

0. Ill 

0. 070 

0.090 

0.072 

1.2 

0.061 

0. 066 

0. 060 

0.051 

0. 062 

0.054 

1.3 

0.045 

0.040 

0.041 

0. 040 

0.042 

0. 035 

1.4 

0.031 

68-4 

0.020 

0.027 

0.011 

0.022 

1.5 

0.023 

60-4 

0.017 

0.021 

48-4 

0.015 

1.6 

0.014 

44-4 

62-4 

0.014 

28-4 

89-4 

1.7 

0. 012 

21-4 

67-4 

0.012 

35-4 

79-4 

1.8 

65-4 

84-5 

18-4 

70-4 

39-4 

50-4 

1.9 

21-4 

28-5 

39-5 

36-4 

63-5 

83-5 

2.0 

14-5 

20-6 

13-6 

62-5 

66-6 

11-5 


The small but significant correlations between 1/a and i and between e and i need 
study, just as the stroi^ interrelation between 1/a and e needs further examination, to 
help understand the nature of meteors. 
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Table 7-11. Distribution of eccentricity. 


e 


Synoptic year 



1961-65 


Observed 

Atmospheric 

Space 

Observed 

Atmospheric 

Space 

0.0 

0.051 

0.169 

0.050 

0.042 

0. 144 

0.021 

0. 1 

0.074 

0.269 

0. 142 

0.052 

0.208 

0.097 

0.2 

0.082 

0.166 

0. 162 

0.059 

0.220 

0. 126 

0.3 

0. 090 

0. 143 

0. 158 

0.067 

0.088 

0. 118 

0.4 

0. 102 

0.082 

0. 142 

0.079 

0. 110 

0. 161 

0.5 

0. 123 

0.078 

0. 135 

0.099 

0. 101 

0. 172 

0.6 

0. 136 

0.060 

0. 120 

0. 122 

0.059 

0. 127 

0.7 

0.139 

0.028 

0. 060 

0. 144 

0. 040 

0. 102 

0.8 

0. 118 

71-4 

0.024 

0. 178 

0. 023 

0. 057 

0.9 

0.076 

33-4 

71-4 

0. 137 

50-4 

0.017 

1. 0 

0.010 

52-5 

13-4 

0. 021 

94-5 

28-4 


Table 7-12. Distribution of perihelion distance. 


q 


Synoptic year 



1961-65 


Observed 

Atmospheric 

Space 

Observed 

Atmospheric 

Space 

0.0 

0. 040 

11-4 

47-4 

0.070 

43-4 

0.014 

0. 1 

0.073 

41-4 

0. 017 

0.113 

0 . on 

0.035 

0.2 

0.067 

77-4 

0. 026 

0. 087 

0. 013 

0.038 

0. 3 

0.070 

0.012 

0.035 

0. 079 

0. 015 

0. 045 

0.4 

0. 070 

0.015 

0. 045 

0.080 

0.024 

0.064 

0. 5 

0.071 

0. 025 

0.055 

0. 070 

0. 028 

0. 079 

0. 6 

0. 086 

0.069 

0. 102 

0.071 

0.084 

0. 114 

0. 7 

0.092 

0.096 

0. 137 

0. 071 

0. 107 

0. 139 

0. 8 

0.129 

0. 187 

0.214 

0.094 

0. 128 

0. 169 

0. 9 

0.245 

0. 464 

0.333 

0.228 

0.496 

0.289 

1. 0 

0.057 

0. 119 

0.032 

0. 037 

0. 090 

0.015 
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Table 7-13. Distribution of aphelion distance. 


q' 

Synoptic year 



1961-65 


Observed 

Atmospheric 

Space 

Observed 

Atmospheric 

Space 

1.0 

0.015 

0.020 

47-4 

0.016 

0.024 

54-4 

1.5 

0.326 

0. 557 

0.434 

0. 263 

0. 501 

0.355 

2.0 

0. 141 

0. 176 

0. 189 

0. 135 

0. 190 

0. 200 

2. 5 

0.104 

0. 083 

0. 106 

0. 105 

0.071 

0. 114 

3.0 

0. 079 

0.048 

0. 073 

0. 088 

0.076 

0.089 

3. 5 

0.066 

0.033 

0.054 

0.071 

0.054 

0.083 

4.0 

0.053 

0.023 

0. 039 

0.053 

0.026 

0. 037 

4. 5 

0. 042 

0.020 

0. 030 

0.043 

0.017 

0. 031 

5.0 

0.030 

0,016 

0.024 

0. 030 

73-4 

0.013 

6.0 

0. 038 

97-4 

0.017 

0.045 

0.011 

0. 023 

7.0 

0. 024 

61-4 

0.011 

0. 028 

54-4 

0.012 

8. 0 

0.016 

14-4 

48-4 

0.020 

55-4 

0.013 

9. 0 

99-4 

10-4 

41-4 

0.012 

22-4 

37-4 

10.0 

81-4 

56-5 

18-4 

0.011 

25-4 

29-4 


0.047 

57-4 

93-4 

0.078 

75-4 

0.019 
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Table 7-14. Distribution of inclination. 


i 


Synoptic year 



1961-65 


(degrees) 

Observed 

Atmospheric 

Space 

Observed 

Atmospheric 

Space 

0 

0. 026 

0. 163 

78-4 

0. 024 

0. 165 

77-4 

2 

0.025 

0. 136 

0.026 

0. 027 

0. 091 

0.021 

4 

0. 024 

0.138 

0.054 

0. 027 

0. 124 

0.040 

6 

0.022 

0. 118 

0.068 

0.023 

0. 121 

0. 062 

8 

10 

0. 020 

0.068 

0.057 

0. 021 

0.066 

0. 037 

15 

0. 048 

0. 163 

0.188 

0. 045 

0. 149 

0. 159 

20 

0.048 

0. 087 

0. 180 

0.046 

0. 093 

0. 150 

25 

0.052 

0.053 

0. 141 

0.048 

0.066 

0. 133 

30 

0. 066 

0. 031 

0.102 

0. 047 

0.044 

0. 116 

35 

0. 060 

0.018 

0. 065 

0. 048 

0. 030 

0.088 

40 

0.061 

99-4 

0.040 

0. 056 

0.020 

0. 064 

50 

0. 123 

91-4 

0.041 

0. Ill 

0.017 

0.066 

60 

0. 124 

37-4 

0.019 

0. 118 

81-4 

0. 035 

70 

0. no 

15-4 

78-4 

0. 108 

35-4 

0.015 

80 

0. 064 

42-5 

23-4 

0.073 

12-4 

53-4 

90 

0.018 

71-6 

37-5 

0.024 

24-5 

10-4 

100 

0.011 

30-6 

16-5 

91-4 

58-6 

33-5 

no 

0.011 

19-6 

12-5 

0.014 

63-6 

31-5 

120 

0.015 

20-6 

12-5 

0.018 

60-6 

27-5 

130 

0.017 

15-6 

84-6 

0.019 

53-6 

16-5 

140 

0.016 

11-6 

64-6 

0. 021 

51-6 

15-5 

150 

0. 019 

16-6 

51-6 

0. 024 

43-6 

76-6 

160 

0.016 

90-7 

20-6 

0. 022 

34-6 

56-6 

170 

0.011 

54-7 

96-7 

0.015 

23-6 

23-6 

180 

58-4 

29-8 

13-7 

99-4 

19-6 

95-7 
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Table 7-15. Distribution of argument of perihelion. 


CO 

(degrees) 

Synoptic year 



1961-65 


Observed 

Atmospheric 

Space 

Observed 

Atmospheric Space 

0 

0.047 

0.053 

0.022 

0.038 

0. 042 

0.020 

20 

0.052 

0.052 

0. 044 

0.050 

0.016 

0. 033 

40 

0. 048 

0.061 

0. 060 

0. 046 

0. 042 

0.041 

60 

0.047 

0. 038 

0.054 

0.044 

0. 049 

0.061 

80 

0. 051 

0.057 

0. 065 

0.045 

0. 092 

0. 071 

100 

0.056 

0. 044 

0. 064 

0.049 

0. 067 

0. 060 

120 

0.073 

0. 072 

0. 095 

0.061 

0.072 

0. 078 

140 

0. 071 

0.075 

0.084 

0.072 

0. 076 

0.075 

160 

0.067 

0.087 

0. 032 

0. 065 

0. 062 

0.035 

180 

0.070 

0. 078 

0. 037 

0. 072 

0. 093 

0. 039 

200 

0.062 

0. 068 

0. 072 

0.069 

0. 102 

0. 066 

220 

0.056 

0.051 

0.083 

0.058 

0.039 

0.073 

240 

0.046 

0. 048 

0. 071 

0. 049 

0.028 

0.071 

260 

0. 047 

0.032 

0.063 

0.050 

0. 031 

0. 069 

280 

0.048 

0. 032 

0.053 

0.054 

0.037 

0.066 

300 

0.050 

0. 026 

0.036 

0.061 

0.030 

0.055 

320 

0. 061 

0. 070 

0. 045 

0.073 

0. 051 

0.062 

340 

0.048 

0.055 

0. 021 

0.045 

0.071 

0. 026 


360 
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Table 7-16. Joint distribution of 1/a and e from the synoptic-year space sample 


l/a\^ 

- 1.0 

9 

- 0.8 
-0.7 
- 0.6 
-0.5 
-0.4 
-0.3 
- 0.2 
- 0.1 
0 . 0 
0 . 1 
0.2 
0.3 
0.4 
0. 5 
0.6 
0.7 
0.8 

0. 9 
1.0 

1 . 1 
1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 
1.9 
2.0 

Sum 


0. i 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 9 1. 0 Sum 

0 0 

0 0 

0 0 

0 0 

2 2 

0 0 

1 1 

1 1 

7 7 

16 16 

52 52 










127 

26 

153 







0 

361 

111 

6 

478 






23 

948 

302 

30 

4 

1306 





2 

1000 

638 

125 

27 

4 

1795 




21 

835 

733 

240 

62 

17 

4 

1912 



3 

847 

822 

245 

120 

56 

16 

5 

2113 


2 

838 

974 

338 

152 

74 

49 

18 

4 

2450 

3 

856 

820 

568 

260 

117 

72 

22 

19 

6 

2742 

548 

964 

635 

329 

183 

101 

67 

33 

9 

5 

2873 

520 

908 

473 

142 

120 

101 

84 

34 

19 

4 

2405 

5 

288 

523 

221 

95 

62 

52 

31 

15 

3 

1295 


44 

201 

240 

197 

75 

71 

41 

15 

4 

889 



0 

73 

159 

95 

41 

34 

18 

6 

425 




0 

56 

164 

70 

61 

22 

5 

378 






41 

65 

16 

12 

1 

135 






0 

59 

62 

17 

6 

144 








12 

21 

6 

39 









5 

3 

8 










0 

0 

1076 

3061 

3493 

3414 

3066 

2909 

2601 

1301 

518 

153 29 

21621 
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Table 7-17. Joint distribution of l/a and e from the 1961^65 space sample 


1/a 



e 


0. 1 0.2 0.3 0.4 0. 5 0.6 0.7 0. 8 0. 9 1. 0 Sum 


- 1.0 

-0.9 

- 0.8 

-0.7 

- 0.6 

-0.5 

-0.4 

-0.3 

- 0.2 

- 0.1 

0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0. 9 

1.0 

1. 1 

1.2 

1.3 

1.4 

1.5 

1 . 6 

1.7 

1.8 
1.9 
2.0 

Sum 


211 

157 

0 


368 


2 2 

2 2 

1 1 

4 4 

1 1 

2 2 

2 2 

2 2 

6 6 

30 30 









0 

55 

55 







0 

225 

44 

269 







434 

142 

20 

596 





0 

440 

350 

81 

12 

882 





1000 

590 

208 

56 

15 

1868 



0 

663 

719 

340 

124 

48 

18 

1912 


9 

521 

643 

307 

170 

93 

44 

16 

1804 

1 

474 

613 

377 

287 

131 

66 

49 

13 

2010 

519 

723 

330 

229 

204 

142 

69 

45 

10 

2272 

652 

396 

297 

359 

113 

92 

62 

35 

11 

2227 

294 

207 

96 

275 

66 

85 

66 

35 

12 

1292 

268 

269 

102 

126 

79 

34 

46 

29 

9 

960 

6 

168 

106 

77 

84 

81 

77 

22 

11 

631 


1 

40 

70 

143 

46 

54 

33 

8 

396 



0 

67 

42 

40 

67 

39 

12 

266 




0 

42 

48 

26 

37 

5 

159 






26 

53 

53 

9 

141 






1 

28 

44 

17 

90 








6 

9 

15 









2 

2 

1739 

2247 

2105 

2886 

3086 

2266 

1822 

1024 

306 50 

17900 
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Table 7-19. Joint distribution of 1/a and i from the 1961—65 space sample. 


cocoo><Dh-co»Ho>CM»-Hc£>CNi^o>eooir-io^^ 

c<ic^^eq<MN>-H^ 


ooooooooooooooooooooo 


ooooooooooooooooooooo 


o o oooooooooooooooooooooo 


oo ooooooooooooooooooo 


O OOOOOOOOO-HOOOOOOO 


oooooooooooooooooooo 


oo ooooooooooooooooooooo 


o OOOOOOOOOOOOOl— 10^0 


O O OOOOOOOOOOOOOr-HOo-HOO 


O O O O O O O O i-H t-H 1— li— (Oi— ICOOCJ O CS] I— ( ' — i O 


OOOoCOTt<COCOil5lOOOOOOO-HOOO^TtCO<No 


oo oooiomooiH(Mi-Hoc<iooiooocoocO'^or-coo 


OOOOOC0Ol0NinOOTt<OO00O(NOOOr^(X)OC0O 


O O OOC4l0c<|OCvIC0lrtT}<0-H00rt<-H(MTiH00C0l0Tt<O(M 

OTflt^OOOOOOcOCQOOJOOTj^Ti^CO O 


O O 0<N or^rJHCOiOrfC^OOOOtOcOOOCDCQOcOTj<c<IC^O 

T^C<IOC^OOOiOr-OOmOOCC»COOO-HT:J<C^ 
00 (M(NMCOC^(NCOO^»H 


ooocoo o Tfincoor-r-oi/5(:Dor-oooior^tr3or-oc<i 

»-HTi<c<icacot^oooior>-i^if 5 >-Ht^oor^(Mo 
OCQCO'^TfOlOOCOC^— I 


o oop-*HTtcacooo«Docaot^OLOOOcaooo> 

I-H «— •COf^lrtOlOOOOOOlrtO'^O |~H t-H 

ocqc<imoocDr-oocq(MrH 


ooccior-cqocaooocaTt^— loorjHococoor^o 
cocor^oo®oiooooQO(Nio 

oiOTfCOCOlOC^C^Oi-H 


oooor^oiOTfcodoooc<icO’^tfDtor'-ooooo<McoTfi/5or-oooo ^ 

>— «000000000000000000000>-H.-H»-5i-H»-H»-HO.-Hcq CZl 

I I I I I I I I I I 


87 


3320 6114 4912 3008 1305 695 289 105 21 7 6 5 3 3 2 1 0 0 19798 


o 

00 


<o(NOQin< 0 C<ir^ 525 S £3 ^ 

^Hc^ino>cJior-Qoo>;^w ^ 

oococoiqco<n^^^* 5 ” ^ 

Cd <M <N C<l C<I ^ 


Si 

a 


o 

o 

Qi 

CQ 

U 

d 

0 
>> 

1 

a 

•l-^ 

I 


rS 


O 

r- 


o 

CD 


o 

ID 


O 




o 

(N 


O 

o 


ooooooooooo o 

OOOOOOOOQOO o 
ooooooooooo o 
ooooooooooo -H 
ooooooooooo ^ 
ooooooooooo ^ 
00000000 »-H 00 cvl 
0000000 --HOOO (N 
OOOOOOOO— lOO CO 


o 

u 

CfH 


o 

o 


o 

oo 


OO-— ^0'~''~•'^'~^’~^'^0 CD 

'-<C<J'^CDLOlDiOmrt<C40 00 

CO 


Q) 


o 


(Noo-^ooor-<Nr^r^co' 

04 ""I C 4 


00 

04 


s 

•r^ 

I 

u 

to 

•rH 

T 3 

I 

o 

>-Z> 


o 

(N 

I 

r-- 

<D 

f-H 

H 


o 

CD 


o 

lO 


o 


o 


o 

04 


>«^co'-Hioior^r-r-o 4 '^»— » 

»-ICO'^lOTfTj<CO<M*^ 


COlDOOr^O*^COCOCOCDC<> 

-H^CDO 04 O 00 r-lD-H 


CDC 40 tOCO’~*COOOOOO'^ 
T:t^CD'^0000CDCOC4r-'(N 
•-H 04 04 D4 <M <M 


T^^OOCDcDCO'-HCDt^COlDO 

CDCDOLOOoO'—trt^'— ICO 
rHri<CDCDlOlO’^ 04 '~‘ 


C0OO00-«t04lD04C0OC0 

O'— •ocooO'^cDOr^'-H 
coooooor^r-o 4 


04 '^r^CO'^OOCDOcDC 4 

O-HO'« 4 <CDiDO 00 Tf 

.-•r-iDinco’cf'co^ 


lO 

CO 


r- 

CD 


r- 

CO 

r^ 


r- 

o 


CO 

04 

o 

CD 


o 

r^^ 


t-H 04 C 0 *'^L 0 cDC^ 00 CJ 5 O'— « 


88 


Table 7-21. Joint distribution of e and i from the 1961 — 65 space sample. 
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8. CONCLUSIONS 


8. 1 Stream Distributions 


We have developed a computer technique for determining the two parameters A 
and CT of the D-distribution of meteor orbits in a stream, and we have applied this to 
the 256 streams found earlier (Southworth and Sekanina, 1973) in the synoptic-year 
data. Moreover, we have relaxed the screening restrictions for high-inclination and 
retrograde streams, discovering an additional 20 streams. We have found that about 
16% of all meteors in the sample belong to the detected streams. 

8. 2 Associations of Streams with Adonis 


Four synoptic-year streams match the orbit of the minor planet Adonis rather 
closely. We undertook to see whether these streams could be understood as meteors 
ejected from Adonis at some time within the last 12000 years, taking into account 
reasonable circumstances of ejection, secular perturbations on the meteors and Adonis, 
and radiation pressure. However, we do not find persuasive evidence for such a 
relationship. 

8. 3 Meteor Heights 

We have improved on our earlier investigation (1973) of meteor heights, revising 
the punched reduction output to a more homogeneous basis of diffusion heights. None- 
theless, we still do not find Ceplecha's "discrete levels" of meteor height, nor any 
important relationships between height and orbital characteristics other than velocity. 

8. 4 Fragmentation 

The spread of fragments of a meteoroid along the trajectory has been computed 
for meteors observed in 1961—65 and in the synoptic year, by using the measured 
number of Fresnel extrema and taking into account diffusion and the limited dynamic 
range of the radar receivers. These spreads will, we hope, constitute some of the 
data for a badly needed physical theory of fragmentation. 
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Our previous conclusion (Southworth and Sekanina, 1973) that fragmentation does 
not constitute a significant selection effect is insecurely confirmed. 

No important relationship was found between fragment spread and orbital elements, 
but a puzzling relationship with velocity has appeared. 

8. 5 Selection Effects 


The statistics Of the length of the ionized column have been revised, taking into 
account the projected extent of the station network on individual trajectories. This 
has yielded new parameters for the relation between column length and radiant zenith 
distance, important both for a physical theory of fragmentation and for a revision of 
our corrections fot selection effects. 

8. 6 Orbital Statistics 

We used the revised corrections for selection effects to update our 1973 orbital 
statistics from the synoptic year and to compute parallel statistics for meteors observed 
in 1961-65. 
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